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and Bruce E. Tinling 
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SUMMARY 


Solar pressure forces can provide restoring torques that aline a 
specified axis of a suitably shaped spacecraft toward the sun-. Implementing 
a control system that employs these torques requires a means of damping oscil- 
latory motion of the solar pointing axis. This paper examines a passive tech- 
nique of damping by dissipating the energy of relative motion between a pair 
of connected bodies. Two systems are studied in which the relative motion 
occurs about a single hinge. Coupling between the attitude motions is caused 
by unequal rates of change of solar pressure torque with each independent 
attitude angle. The performance is compared with that of a third system that 
uses a two-degrees -of -freedom hinge and has no coupling between the attitude 
motions . 


INTRODUCTION 


Solar radiation exerts a pressure which can cause torques large enough to 
be an important consideration in the attitude control of spacecraft, partic- 
ularly when the spacecraft mission is of long duration. In most instances 
solar torques cause undesirable disturbances that must be countered by the 
control system. However, for spacecraft required to point continuously toward 
the sun, solar torques can be an asset if the spacecraft is shaped so that the 
desired attitude orientation is also one of stable equilibrium with respect to 
the solar torques. Both passive and semipassive control systems that rely on 
this principle have been proposed (refs. 1, 2, and 3). 

Attitude control systems that exploit solar pressure torques differ 
primarily in the technique used to damp the attitude motions. In some, the 
damping may be active, as in the control system for the Mariner spacecraft 
(ref. 4). In others, the damping may be passive, requiring no sensors, actu- 
ators, or additional energy sources. An example of a system of this latter 
type is described in reference 3 where a structural element provides the phase 
shift between the solar restoring torque and attitude motion necessary to damp 
the attitude motions. 

An alternate passive technique for solar oriented attitude control 
systems is proposed in reference 5 • This scheme aims to remove any unwanted 
mechanical energy, resulting from external disturbances, by connecting to the 
spacecraft, through dampers, one or more auxiliary bodies. If the 



configuration of these auxiliary bodies is correct, any attitude motion will 
cause relative motion between the spacecraft and auxiliary bodies, and the 
energy contained in this relative motion will be absorbed by the dampers. 
Mechanical simplicity suggests the desirability of using a single auxiliary 
body with a single degree of freedom relative to the spacecraft. To achieve 
this simplicity it is necessary to promote coupling between the degrees of 
freedom of the entire system so that any attitude error will cause relative 
motion between the two bodies. Two techniques of coupling the modes of motion 
are possible. One is to distribute the mass asymmetrically in either or both 
bodies. The performance of a solar oriented spacecraft utilizing this tech- 
nique has been examined and reported in reference 5. The other technique is 
to create unequal rates of change of solar torque with each independent atti- 
tude angle by means of suitably shaped illuminated surfaces of either or both 
bodies. Any attitude deviation of either body then produces a solar torque 
that causes relative motion between the two bodies. This technique is demon- 
strated in this report by an analysis and comparison of two possible single- 
hinge configurations, both exhibiting ,T torque coupling." 

The performance of this type of solar oriented spacecraft probably can 
be improved by increasing the number of bodies or the number of degrees of 
freedom between the bodies. The price for the improved performance is addi- 
tional mechanical complication. To illustrate the magnitude of the perfor- 
mance gains obtainable in this manner, a third system consisting of a single 
auxiliary body having two degrees of freedom relative to the main spacecraft 
body is analyzed. This system is made as simple as possible by arranging for 
both bodies to have symmetrical mass distributions and solar reflecting sur- 
faces relative to the sun line and two hinges orthogonal to each other and to 
the sun line. The attitude motions of such a system must be uncoupled. 

The study is restricted to heliocentric orbits wherein the natural 
frequency of oscillation of the spacecraft is large compared with the orbital 
frequency. This implies that the torques due to the gravity gradient of the 
sun and the gyroscopic effect of the curvature of the orbital path are negli- 
gible compared to those derived from the solar pressure. It follows that 
solar stabilization can only be achieved about two axes normal to the sun line. 
It is assumed that the spacecraft is stabilized to prevent rotation about the 
sun line by other means . 


NOTATION 


Coordinate Systems 

cT.(j - 1,2,3) unit vectors of the vehicle hinge coordinate system defining 
J a reference frame fixed in the main satellite body and 

derived from the y.(j = 1,2,3) reference frame by rotation 

about v 3 through the hinge angle 7 (The hinge axis is 

) 
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^ j ( j — 1,2,3) 


Oj(J = 1,2,3) 


▼j(J = 1,2,3) 


unit vectors of the damper coordinate system along principal 
axes of the damper body (in the undisturbed steady state 

d. = v (j = 1,2,3) 

J J 

unit vectors of the orbital coordinate system defining a 
right-hand orthogonal reference frame fixed in the 
satellite orbit such that o 3 is directed toward the sun 
(Orientation of the reference frame about o 3 is 
arbitrary. ) 

unit vectors of the vehicle coordinate system defining the 
torque decoupling and principal axes of_the main satellite 

body (in the undisturbed steady state v. = o.(j = 1,2,3).) 

J J 


Symbols 


B 


C 


dimensionless damping constan 


dimensionless spring constant 



D viscous damping constant of angular rate damper located on 

c 2 axis 

I (j = 1,2,3) principal moments of inertia of main satellite body 

V j 


K 


J, K 

di a 2 


K , K 
v x > v 2 


L 


spring constant of spring restraining rotation of damper 
body relative to main body about cT 2 axis 

solar spring rates of damper body about d x and d 2 axes, 
respectively 

solar spring rates of main body about v ± and v 2 axes, 
respectively 


undamped natural frequency parameter 


(l 4lV\ 

V.KV! la J 


s 

s 1 



Laplace transform variable expressed in dimensionless time 
Laplace transform variable expressed in real time 
real time 

external torque acting on damper about its hinge 


3 



T. (j = 1*2,3) external torques, other than solar torques, acting on 

complete satellite 

— vector rotation of damper body with respect to inertial 

space 

■jr vector rotation of main body with respect to inertial 

p space 


7 hinge angle (See definition of c.) 

J 

6 imaginary part of a root of the characteristic equation 


Q ^ small rotation angle of damper about c 2 axis 

(■ vector rotation of damper body with respect to main body 


a 


T 


real part of a root of the characteristic equation 


dimensionless time, 



cp, 0, small rotation angles about v x , v 2 , and v 3 relating the 

v . frame to the o . frame 
3 3 


c, h 


d 

d.(j = 1,2,3) 

S 


v 

Vj(d = i> 2 >3) 


Subscripts 


quantity measured with respect to the 
respectively 

damper body 

quantity measured with respect to the 
solar torques 
main body 

quantity measured with respect to the 


c or h frame, 


d . frame 
3 


v . frame 
3 


Double Subscripts 

vc, etc. variable used in describing relationships between the 

v. frame and c^ frame, etc. For example 

is direction cosine matrix relating the v frame 
relative to the c frame 
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SCOPE OF ANALYSIS 


The analytical technique adopted is aimed at establishing the broad 
feasibility of this type of system and at comparing the performance of three 
configurations. Unlike the performance of many spacecraft control systems, 
the steady-state performance of those considered here is not dependent on any 
presently well-defined external disturbances. Apart from the meteorite prob- 
lem, possibly the most significant limitation on steady-state pointing accu- 
racy is the inhomogeneity of surface solar reflection properties due to 
nonuniform aging, and inaccuracies in the torque-free, steady-state alinement 
of the two bodies due to manufacturing errors. Since these considerations are 
not strongly related to the problem of selecting fundamental parameters, the 
question of steady-state performance is not treated in this report. 

The remaining significant control system performance criterion is the 
damping time, which dictates how long it takes the system to reach the steady 
state after injection or after some subsequent disturbance. Consideration of 
the damping time goes far toward establishing the bounds on many of the impor- 
tant design parameters for this type of system. As in reference 5, the analy- 
sis is based on the linear autonomous equations approximating the system 
behavior. The damping criterion is the damping of the least damped mode of 
motion, which, in this case, is characterized by the maximum real part of the 
roots of the characteristic polynomial. Minimizing this quantity is then 
equivalent to insuring that any arbitrary attitude disturbance decays at the 
maximum rate. For the coupled configurations, this minimization is accom- 
plished by the steepest ascent computational techniques described in refer- 
ence 5* An analytical solution for the minimization exists for the two -hinge 
uncoupled configuration. (See appendix B.) 

In principle, it is possible to find the best damped two-body, one- or 
two -hinge configuration without restricting the mass distribution, body shape, 
or hinge orientation. In practice, the computational problems and, in partic- 
ular, the slow rate of convergence to a solution, make such a program too 
ambitious, at least for the present. The best that can be attempted at this 
time is to choose certain of the parameters so as to promote some symmetry, 
thereby reducing the computations to a more manageable level. In particular, 
it is assumed throughout that the body geometry and mass distribution permit 
the existence of "torque decoupling" axes as described in appendix A and that 
these are coincident with the principal axes of inertia. Furthermore, the two 
bodies are assumed to be connected so that their centers of mass coincide and 
their relative orientation is assumed to be such that, in torque-free equilib- 
rium, their "torque decoupling" axes and therefore, also their principal axes 
of Inertia, coincide. These restrictions imply the existence of two orthog- 
onal planes relative to which the mass distribution and surface geometry of 
both bodies are symmetrical. When the intersection of these two planes is 
alined with the sun line, no torques are acting on either body. The hinge or 
hinges are assumed to be perpendicular to the intersection of the two planes. 
One of the .chief implications of these assumptions is that inertial coupling 
due to asymmetrical mass distribution is excluded. This type of coupling has 
already been investigated in reference 5 ■ 
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It is pointed out in reference 5 that the damping of passive solar 
oriented systems, expressed in terms of real time, can always he improved by 
increasing the natural frequency. This, in turn, might imply unrealistic 
ratios of the solar spring constant to the moments of inertia of the satellite 
bodies. To make the solutions physically realizable and yet avoid specific 
practical design constraints, the damping is normalized by expressing it in 
terms of a unit of ti me propor tional to the undamped frequency in roll of the 
main satellite body n /k Vi /I v - 1 , where is the solar spring rate and I Vi 

is the moment of inertia. In addition, the point of view is adopted that the 
main sate llite body will be constructed to have the largest practical value of 
and, therefore, with all other parameters sized accordingly, the 

maximum damping, and that the remaining corresponding parameters n/K ( j 1 /I ( j 1 

n/|Kv 2 1/Iv 2 will be no greater than this* The argument in favor 

of this i s that if any of the above parameters could be made larger than 
n /Ky 1 /I Vi ^ then the same constructional technique could be used to increase 

n/K Yi /I Vi to the same value. This idea can be expressed in the form of the 
inequality 



The linear equations of motion of a class of configurations sufficiently 
general to include as special cases all those configurations of specific 
interest are derived in appendix A 


ANALYSIS OF SPECIAL SPACECRAFT CONFIGURATIONS 


Three separate spacecraft arrangements were studied -which are classified 
according to the variation of solar pressure torques with attitude. Two 
types, described as the wedge-wedge and cone-wedge configurations, have 
unequal rates of change of solar pressure torque with each independent atti- 
tude motion. These two types, therefore, can be damped when a single degree 
of freedom of the relative motion exists between the coupled bodies. The 
third type, described as the cone-cone configuration, is axially symmetric 
about the solar pointing axis in both mass distribution and external shape. 
The cone-cone configuration, therefore, requires a two-degrees -of -freedom 
hinge. These three configurations will be discussed separately in the 
following sections. 


Wedge -Wedge 

This configuration is probably the simplest conceivable, two -body single - 
hinge, sun-stabilized system. Its essential features are illustrated in fig- 
ure 1. Both the main satellite body and damper body are shown as wedges to 
indicate that they have solar restoring torques about one axis only. In 
torque -free equilibrium, these "solar spring" axes are perpendicular to each 
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Figure 1.- Wedge-wedge configuration. 


other and, together with the hinge line, 
are perpendicular to the sun line, or yaw 
axis. The roll and pitch principal iner- 
tias of the main satellite body are 
equal. These special characteristics 
can he summarized as follows: 



K .. = 0 

<ii 

I = I 

vi v 2 J 


( 1 ) 


The characteristic polynomial of this 
system may he obtained directly from the 
linear dynamical equations derived in 
appendix A. Using conditions (l), the 
polynomial, in dimensionless form, is 
the determinant 



(s 2 +l)sin 7 


0 




«a£vA 

Kv 1 Id 2 / 



-s 2 cos 7 


- ^ — s^ sin 7 



Kd 2 IvA 

+ K Vl IdJ C ° S 


7 


(Bs+C) A 

± Vl 


= 0 


( 2 ) 


where s is a dimensionless Laplace transform variable and 
C = (K/I c i 1 ) (l Vl /K v ) and B = (d/I^) JK^/T^ are the dimensionless hinge 
spring and damper constants. The factor s 2 = 0 in equation ( 2 ) merely repre- 
sents the fact, mentioned earlier, that motion about the sun line (yaw axis) is 
not stabilized. It is clear from equation ( 2 ) that the hinge axis must be 
skewed to the torque decoupling axes because if sin 7 = 0 or cos 7=0, an 
undamped mode of motion is present. Thus, a system with damping about both 
pitch and roll can be achieved only if 


7 4 -f n = 0,1,2 (3) 

In addition, the ratio of the constant coefficient to the coefficient of the 
highest power (s 6 ) of the characteristic polynomial is 


7 


[ 2 


Iv! V K v x , 


x d 2 A&i 


^Vx Vjvj 


+ 1 Jcos^ 7 + 


I(3 -i 

*^ v i V^vx 


(4) 


+ 1 Jsin^ 7 


and since, for stability, this quantity must be positive, it follows that 
Cfe^/Kvx) (idxAvx) mus t t>e positive. Since C = (K/K Vl ) (ivxAdx) it follows, 
therefore, that K/K & 2 must be positive for stability. This requires that 
K<} 2 be positive since K must also be positive for a mechanical spring. In 
other words, the solar torques on the damper must be in the stable sense. The 
above derived conditions necessary for stability are illustrated in figure 1 
80( _ by the position of the hinge line and 

the shape of the damper. 
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Figure 2.- Optimum damping and associated 
values of free parameters for the wedge - 
wedge configuration. 


The results of computations by the 
method of steepest ascent (ref. 5) to 
find the best damped wedge -wedge con- 
figurations are shown in figure 2. The 
damper inertia parameter I^/l^ has 
been chosen as the independent variable, 
and the other free parameters shown are 
required to achieve optimum damping for 
a particular value of 

damping parameter 


maio 


idi/i 


vi* 


Vl 


The 

repre - 


sents the real part of the root of the 
least damped mode normalized by the 
natural frequency of the main body 
alone . 

As noted in the Introduction, the 
performance of an alternative scheme 
which relies on asymmetric mass distri- 
bution to promote the necessary coup- 
ling for a single -axis damper has been 
reported in reference 5* A comparison 
of results indicates that a comparable 
spacecraft of the wedge -wedge type 
would be better damped by a factor of 
3 or more. 


Cone -Wedge 

Some improvement in the performance of the torque -coupled, single-hinge 
systems might be expected if the main body produced restoring torques about 
all axes normal to the sun line. A special configuration of this type was 
studied in which the main body was assumed to be axially symmetric in mass 
distribution as well as in external shape. As in the wedge -wedge arrangement. 
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the hinge axis was assumed to be in a plane normal to the axis of symmetry 
The special conditions for the cone-wedge can be summarized as follows. 


K = K 
Vi v 2 

K, = O' 
o-i 

I = I 
Vi v 2 


( 5 ) 


Using conditions ( 5 ) , the determinant of the characteristic polynomial of the 
dimensionless linearized equations of motion becomes: 


(£r +1 > 


■+1 


E vi 


s^ sin 7 


(s 2 + l) 


sin 7 


(wf + + (& i&) 

-(s 2 +l)cos 7 




L n. 


Ek 

Z Vx 




-bos y 


ia, 

(C+Bs) — 




= 0 


(6) 


As in the case of the wedge-wedge system, an undamped mode of motion exists 
unless the hinge axis is skewed to the torque decoupling axes, so that the 
condition expressed by ( 3 ) must be satisfied. The ratio of the constant coef- 
ficient to the coefficient of the highest power (s 6 ) of the characteristic 
polynomial is , in this case. 


-C 1 + 


K d 2 \ I dx K ds 

i? — cos^ 7 
h Vl 


K- 


•Vi./ - L v 1 


I d 2 /i 




Icos 2 7 


Idi ( lds A A • 

- = T + 1 PI 


L Vx \ Vx 


L Vx V-Vx 




is in 2 7 


which must be positive for stability. This implies that 


Kc 

K- 


Vx 


cos 2 7 


(&) 


+ c 


( 7 ) 


(8) 


Since C is positive for a real spring in which torque opposes rotation, the 
ratio of the solar "springs" and K Vl must always be more positive than 

some negative number. Since the manner in which the Laplace variable was 
normalized implies that K^ 1 is always positive, it follows that the crite- 
rion can be satisfied when the solar torques on the damper body are destabi- 
lizing as well as when they are stabilizing. This possibility is illustrated 
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Figure 3 - Cone -wedge configuration. 


80 


by the cone -wedge arrangement in figure 3 
where the wedge , representing the damper 
body, is shown open toward the sun. 

The results of the computations to 
find the best damped cone -wedge systems 
are shown in figure k. The damper body, 
for all damper inertia ratios, is shown 
to be destabilized by the solar pressure 
torques (K^/K^ negative). This system 
should be better damped than the wedge- 
wedge because the torques on the coupled 
bodies caused by a given attitude error 
will be opposite in sense, thus leading 
to the maximum relative motion and, there- 
fore, to the maximum energy dissipation 
in the damping mechanism. This is borne 
out by the results which show that the 
damping for the cone-wedge system is 
about 20 percent greater than for the 
best system of the wedge -wedge type. 
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Figure - Optimum damping and associated 
values of free parameters for the cone- 
wedge configuration. 


Cone -Cone 


When both connected bodies are 
axially symmetric, the small angle 
motions are not coupled and, therefore, 
two degrees of freedom are required. 
Axial symmetry requires that: 


K 

Vl 


K 

v 2 



I 

■Vi 



I 

v 2 




( 9 ) 



When these conditions are imposed and 
the hinge angle / is chosen to be tt/2 
the linearized equations of motion given 
in appendix A yield the following for the 
determinant of the characteristic 
polynomial: 
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+ 1 S' 



+ 






( 





(s 2 + 1) 




Sun-line 


Figure 5** Cone -cone configuration (two- 
degrees -of -freedom hinge) . 



Figure 6.- Optimum damping and associated 
values of free parameters for the cone- 
cone configuration. 


As would be expected, there are two 
modes of motion which are undamped- One 
of these modes represents rotations about 
an axis normal to the hinge axis and will 
be damped if an additional degree of 
damper freedom is allowed. 

The optimization of the damping of 
the planar oscillations represented by 
the 2x2 determinant of equation (10) can 
be cast in a form analogous to the damp- 
ing of the pitch librations of a gravity 
stabilized satellite such as that ana- 
lyzed in reference 6. Pertinent results 
from reference 6 are used in appendix B 
to derive the parameters for the best 
damping for two-degrees-of -freedom solar - 
stabilized spacecraft. A typical cone- 
cone configuration is sketched in 
figure 5 and formulas for the optimum 
damping and associated parameters are 
given in table I- 

The results for two-degrees -of - 
freedom dampers are given in figure 6. 

The damping exceeds that for the single- 
hinge systems by a factor of 2 or more. 

To achieve this damping, the damper body 
is shown to be destabilized by the solar 
pressure torques. From Avi - 0 to 
0.382, (K dl /l Vl )/(K Vl /l Vl ) = -1, indicat- 
ing that, per unit of inertia, equal and 
opposite torques are generated on each 
body in response to a disturbance. As 
will be noted in the next section, desta- 
bilizing solar pressure torques on the 
damper body are undesirable in some 
instances. The performance was also 
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studied, therefore, for the case when the damper is specified to have no solar 
torques (K^ =■ 0) . For this case the damping is approximately halved. 


DESIGN CONSIDERATIONS 


Most missions for which a passive, solar -stabilized spacecraft is attrac- 
tive are likely to be extremely long and may require an elliptic heliocentric 
orbit. Both of these factors can cause significant departures from those con- 
ditions for which the damping of the attitude motion has been maximized. On 
lengthy missions the properties of the spacecraft surfaces are expected to 
deteriorate, thereby changing the restoring torques. A more pronounced, and 
predictable, effect will be caused by elliptic orbits. These orbits cause off- 
optimum conditions because of the inverse square variation of solar pressure 
with distance to the sun. The criteria for ' selecting a particular’ design, 
therefore, will probably be heavily weighted toward maintaining adequate damp- 
ing despite departures from the design parameters. 


1000 


100 


t sd=i 


io 





I 2 


Cone -cone (optimum) 

Cone -cone (K<j = 0) 

Wedge-wedge 

Cone -wedge 



Figure 7*- The effect of solar distance on 
damping. 


The effect of variation of distance 
from the sun is shown in figure 7 for the 
configurations that have been studied. 

The results are presented to show the 
effects of departure from the solar dis- 
tance for which the performance is an 
optimum. The ordinate Tg^_^ is the 

time to damp to l/e of the initial 
amplitude normalized by the natural fre- 
quency when SD =■ 1. The curves, there- 
fore, represent the variation of damping 
time with relative distance from the sun. 
Several of the configurations are shown 
to be unstable when the solar distance is 
decreased by as little as 10 percent. 

This instability is traceable to the 
presence of a destabilizing damper body. 
Although such bodies consistently provide 
the best damping for optimum conditions, 
they are more sensitive to off -design con- 
ditions • Instability results when the 
destabilizing torque on the damper over- 
comes the restoring torque on the mechan- 
ical spring that positions the damper 
relative to the main body. The damped 
body then moves to one side until the 
limit is reached, causing an attitude 
error and cessation of damping. 


The instability with decreasing solar distance is avoided if a configura- 
tion is chosen for which the damper body is at least neutrally stable with 
respect to solar pressure torques. For the wedge-wedge configuration, the 
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damper body is required to be stable and the damping does not diminish drasti- 
cally as the solar distance is decreased. The cone-cone arrangement can be 
specified to have no solar pressure torques on the damper. This arrangement 
has somewhat better performance than the wedge-wedge for almost any solar 
distance . 

The choice of a single-degree-of -freedom or two -degrees -of -freedom damper 
would most certainly be based on considerations other than the level of damping 
afforded by each. For either damper, the crucial item in the mechanization of 
a passive solar stabilization scheme is most likely to be the spring -damper 
mechanism which connects the main and damper bodies. When it is realized that 
the whole scheme depends upon pressures, which at Earth distance are on the 
order of only 0.1j-5xl0“ 5 newtons /meter 2 , it becomes apparent that the spring- 
damper must be entirely free of friction, or "stiction" as it has come to be 
called, and must furnish a precise, but small, spring constant. Spring -damper 
devices with these characteristics have been designed and built for passive- 
gravity -oriented spacecraft. (See refs. 7 > 8* and 9 -) Of particular interest 
for application to solar spacecraft are devices that employ a diamagnetic 
bearing . 


CONCLUDING REMARKS 


Several techniques have been studied for passively damping an inter- 
planetary spacecraft that is oriented toward the sun by solar pressure torques. 
All of the methods considered rely on dissipating the energy of the relative 
motion between two connected bodies. The principal difference in the tech- 
niques was whether one or two degrees of freedom of the relative motion was 
permitted. When only a single degree of freedom was permitted, it was neces- 
sary to provide coupling between the attitude motions so that any disturbance 
would cause relative motion between the two bodies. In this study, the coup- 
ling was provided by different rates of change of torque for each independent 
attitude angle- 

Although the damping of the single -hinge spacecraft was always inferior, 
single-hinge arrangements were found that were sufficiently well damped to 
reduce the oscillation amplitude by 60 percent within the time required for 
one cycle of oscillation of an undamped spacecraft. It is possible to provide 
better damping for either arrangement by shaping the auxiliary, or damper, body 
to cause the solar pressure forces to be destabilizing. The damping of such 
arrangements, however, is extremely sensitive to changes in the distance from 
the sun and have an unstable mode when the solar distance is decreased by as 
little as 10 percent . 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif., 94035 , Mar. 28, 1968 
125 -19-03 - 04 - 00 -21 
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APPENDIX A 


DERIVATION OF THE LINEARIZED EQUATIONS OF MOTION 


The linearized equations of motion are developed for a spacecraft config- 
uration sufficiently general to include, as special cases, all the spacecraft 
configurations analyzed in this report. The conditions defining the degree of 
generality of the configuration along with other simplifying assumptions are 
as follows : 

(a) A set of three mutually orthogonal axes (torque decoupling axes) can 
be found for both the main satellite body and damper body, such than an angu- 
lar rotation about any one axis produces no solar torques about either of the 
other two axes . 

(b) The undisturbed steady-state orientation of the main satellite body 
and damper body relative to one another is such that their torque decoupling 
axes coincide. Furthermore, in undisturbed steady state one of these axes 
coincides with the sun line, and angular rotation about this axis produces no 
solar torques about either body. 

(c) The vector solar torque acting on either body in the undisturbed 
steady-state orientation is zero. 

(d) The principal axes of inertia and torque decoupling axes of the main 
satellite body and the damper body coincide. 

(e) The centers of mass of the damper body and main satellite body 
coincide. 

(f) Both the main satellite body and damper body are rigid. 

(g) The orbital angular velocity of the satellite is small compared with 
the circular frequency of oscillation of the satellite. 

The linearized dynamical equations of the satellite system are derived 
most elegantly from Lagrange’s equations of motion. 

For small rotations, the orientation of the main body and damper from 
their initial positions can be described by a vector rotation quantity. Thus, 
the vector rotation p of the main body with respect to its equilibrium 
position is 

P = |vi, v 2 , v 3 | cp 

0 (Al) 

♦ 


Ik 


The rotation £ of the damper body with respect to the main body is 


I “ | c l^ C 2 > °3 


0 

e. 


= |Vl, Vfi, V 3 | My-Q 


where to a first order 


M 


ve 


cos 7 
sin 7 
0 


-sin 7 0 

cos 7 0 

0 1 


Substituting equation (A3) into equation (A 2 ) produces the following. 


I = |vi, v 2 , v 3 


- 0 d sin 7 

Q cos 7 
d ' 


(A2) 


(A3) 


Therefore, the rotation a of the damper with respect to inertial space is 

a = | + P = | Vi, v 2 , v 3 | 


cp - 0 d sin 7 

0 + 0, cos 7 
d ' 

■f 


(A 4 ) 


The kinetic, potential, and dissipative energy expressions may now be written 
as follows 


T = — 

2 






I, a, 2 + I m 2 + I 0 2 + I i 2 
dj j vi r v 2 v 3 


K , a 2 + K cp 2 + K 0 2 + K i 2 + K0? 
d. j vi r v 2 v 3 y d 

J 


(A5) 

(A6) 

(AT) 
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Lagrange's equations for the system are. 


_d_ 

fbT) 

, 8V 

dt 

\§cp / 

+ "T-™ + 

ocp 

d 


+ dv + 

dt 

KdeJ 

+ 80 + 

± 

/oT) 

+ $z + 

dt 

w 


d 


^ + $L- 

dt 

V^d' 

/ 80 d 


= T 


L 2e 


L 3e 


= T. 


(A8) 

(A9) 

(A10) 

(All) 


Substituting for T, V, and A from equations (A5), (A6) , and (A7) into equa- 
tions (A8), (A9) , (A10), and (All) and using equation (AU) for a yields the 
following equations of motion. 


I, + I cp + ( K , + K cp - ( I , § , + K, 0, Is in 7 

<*1 v/ V d x r V did d/ 


+ )S + ( K d2 + K v je + U d2 8 d ♦ K d2 e a j=os r 


= T 


ie 


= T, 


2e 


I, + I H 

d-3 v 3y 


= T, 


3e 


(A12) 

(A13) 

(Al4) 


-^d^ + K d 1 CP ) Sln 7 + ^d/ + K d 2 0 ) C ° S 7 


+ ( ! da cos- 7 + I di sin- y)9 a 


+ (K da cos 2 7 + K dx sin 2 7 fa + K0 d + D0 d = T 4e 


(A15) 


The characteristics determinant of this system is 
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-S. 


(vO'XvO 


in 7 


(v^'XvO 


-(v**0' 

(v’"*0” 5 ? 


(> 


I. +1 
<33 v a 


- ( I a 1 s,2+K d 1 ) sin 7 (v’^O 01 


:OS 7 


(> 


I, cos 2 7+1, sin 2 7 
d s o-i 


+(K, cos 2 7+K, sin 2 7+K 


^'Xdg' 

=) 


= 0 


(A16) 

By standard operations this determinant may he reduced to the following form. 


(’ 


I , +1 
d3 v 3> 


(vv>' 2 * (vO 


I s T +K jsin 7 
Vi V 3 ./ 


0 

7 

(vO'* + (vO 

(v’^v) 5 ” 7 

Vv 7 

(Ds’+K) 


= 0 


(A 17 ) 


Equation (A17) may be made dimensionless by replacing s’ by and 


dividing each row by . The result is 


(Sr + 0 s2 + (Sr 


+ 1 


(s 2 + l)sin 7 


-(£ 


Isin 7 


os 7 


20s 7 


Bs + C) =-*■ 
Iv x 


= 0 


(A18) 


where constant factors have been deleted. 
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APPENDIX B 


OPTIMUM DAMPING OF TWO -DEGREES -OF -FREEDOM UNCOUPLED CONFIGURATION 


It follows from equation (10) that the characteristic equation for the 
motion about the Ox axis of an uncoupled two -degrees -of -freedom configura- 
tion is 


s 4 + s 3 B 


*31 

: vi 


+ 1 + s‘ 


dl 

•Vi 


+ 1 + 1 + L 


+ sB 


X d! 

L + 1 

IVi 


C L 




+ 1 + L 


- 0 


(Bl) 


where 


L - 




and |Lj ^ 1 


Equation (Bl) has the identical form of the characteristic equation 
treated in reference 6. Rather than attempt to interpret the solutions given 
in reference 6 in terms of the particular notation adopted here, one of the 
general results of reference 6 will be used to derive the particular solutions 
of interest. 

It is shown in reference 6 that the best damping for the system under 
consideration occurs when the characteristic equation has repeated pairs of 
complex roots. The characteristic equation therefore can be written in the 
form 


(s + of + iS) 2 (s + CT - ib) 2 = 0 (B2) 

where cr and 5 are real quantities. The expanded form of equation (Bl) is 

s 4 + has 3 + (2S 2 + 6cr 2 )s 2 + ba(a 2 + 5 2 )s + (cr 2 + 5 2 ) 2 = 0 (B3) 

The coefficients of equations (Bl) and (B3) may be equated; thus. 
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4<J = B 

(r t+ 0 

(b4) 


Vp-vi / 

2(6 2 + 3 a 2 ) = C 

0 +ltL 

(B5) 

4 cr(cr 2 + 6 2 ) = B 

(l^+i) 

(B 6) 


\ I v , l / 


(a 2 + S 2 ) 2 = C 

( L ^r +1 ) +L 

(B7) 


The quantities 5, B, C may be eliminated from equations (B4) to (B7) to 
yield the following expression for a in terms of 1^ /i^ and L. 



Dividing equation (B6) by equation (B4) yields 


1 + L 


a 2 + 5 2 


L di 


L vi 


1 + 


E Vl 


(B9) 


and substituting for a from equation (B8) into equation (B9) results in the 
following expression for S 2 


k 1 + L 


- 


v 2 


C vi 


(1 - L Y‘ 


4 1 + 


C vi 


1 + L 


E Vl- 


(BIO) 


A necessary condition that the solution be of the form given by equation (B2) 
is that 
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cr 2 2: 0 
& 2 :> 0 

From equations (b 8) and (BIO) these conditions become 


x d 1 ' N 

1 + L - — ) :> 0 
I Vl 


, W Ida , 

1 *( 1 + L ir-) - — t 1 - L > £ 20 




L V! 


It follows from equation (B8) that for a given value of 1^ Avi the smaller 
the value of L the greater the value of the damping coefficient a. Thus, 
maximum damping for a given value of I^/l^ occurs at the smallest value 
of L for which the following conditions hold: 


1 + L 


V 2 


L V 1 


Zi' 

r Vl 


di 


£ 0 


M 1 + L J — ) - r (l-L) 2 kO 




(Bll) 


l| * l 


The third of the inequalities (Bll) shows that the lowest bound on L is -1. 
From equation (B8) the corresponding maximum damping is 




max 


11 - 


di 

•Vi 


(B12) 


However, the range of validity of this solution is given by the first two of 
inequalities (Bll) which, when L = -1, become 

■'V 


Wl 


£ 1 


dA 

-) -3r+i^o 

•Vl / X V! 


(B13) 
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The second of inequalities (B13) holds if, and only if 


£ ^ ^ and 


■^1 




3 

2 


However, the first of inequalities (BI 3 ) excludes I dl /l Vl ^ (3 +*/5)/2 and 
since physical considerations ensure that "the solution (B12) 

exists if, and only if, 


0 



3 -JT 

2 


(Blk) 


When I ( j 1 /l v - 1 > (3 -n/ 5)/2, the first two of inequalities (Bll) must determine 
the minimum value of L and, therefore, the maximum damping. Since, when 
I dl /l Vl = (3 - >J3 ) / 2 and L = -1 the solution lies on the boundary of the 
region defined by the second of inequalities (Bll), it seems reasonable to 
examine the solution along the boundary for I d /l Vl > (3 -*[5)/2. In this 
case, the second of inequality (Bll) becomes an equality and may be solved for 
L 


L 



(B15) 


Substituting this value of 


L into the first of inequalities (Bll) yields 



which is always greater than zero when I ( ^ i /l Vl >0 so that the value of L 
which satisfies S 2 - 0 satisfies the first of inequalities (Bll) when 
IdiAvi > 0. Therefore, the solution on the boundary 6 2 holds when 

Idi/lva. > (3 “n/ 5)/2 and the value for the maximum darling from equation (b 8) 
is 


max 


L di 


1/2 


(Bl6) 



Graphs of a mg:x against I^/l^ from equations (B12) and (Bl6) are given in 

figure 6. The corresponding values of B and C can be derived from 
equations (b 4) to (B7), 


C 


B = 


ka 

max 



(1 





(B17) 


(B18) 


If the damping body is constrained to have no solar torques acting on it, then 
L = 0. For this condition the first and third inequalities of (Bll) are 
always satisfied, while the second of inequalities (Bll) is satisfied if 






£ 


4 


Therefore, in the region 0 £ I ( ^ i /l Vl £ 4 the maximum damping is given by 



The damping constant B and spring constant C are given by equations (B17) 
and (B18), respectively, with L = 0. When 1^ Avx >4, the solution on the 

S 2 = 0 boundary holds and J > 0. 


A summary of all pertinent formulas is given in table I 
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TABLE I.- FORMULAS FOR UNCOUPLED SYSTEMS 







